Capillary Rise of Liquids in Nanopores 
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We present measurements on the spontaneous imbibition (capillary rise) of water, a linear hy- 
drocarbon (n-Ci6H34) and a liquid crystal (80CB) into the pore space of monolithic, nanoporous 
Vycor glass (mean pore radius 5 nm). Measurements on the mass uptake of the porous hosts as 
a function of time, m{t), are in good agreement with the Lucas- Washburn v^-prediction typical of 
imbibition of liquids into porous hosts. The relative capillary rise velocities scale as expected from 
the bulk fluid parameters. 



I. INTRODUCTION 

Spontaneous imbibition, the capillary rise of a wetting 
liquid in a tube or more generally spoken a porous host, 
plays a crucial role in such different areas as oil recovery, 
printing, cooking and fluid transport in living organisms 
[ij . It is governed by a balance of capillary action versus 
viscous drag and gravity forces. Therefore, the details of 
the imbibition process, e.g. the morphology and speed 
of the advancing imbibition front depend sensitively 
on fluid parameters (viscosities 77 of displacing and 
displaced fluid and interfacial or surface tensions cr) and 
on the structure of the porous host. 

Here, we shall present experiments on spontaneous 
imbibition of liquids into nanoporous, monolithic silica 
glass. We will demonstrate that such experiments allow 
one to explore flow properties of liquids in restricted 
geometry. A topic which is of importance for lubrication, 
transport of liquids through biomembranes and folding 
of proteins - not to mention the growing, technological 
and scientific interest concerning the physics of liquids in 
nanodevices, e.g. carbon nanotubes [3) S S H, H, 0, Hi- 

As porous host we have chosen "thirsty" Vycor glass, 
a nanoporous silica glass which is permeated by a 3D 
network of interconnected, randomly oriented pores 
In order to play with the complexity of the imbibing fluids 
we have chosen three different liquids: water, a chain- 
like, linear hydrocarbon (n-Ci6H34) and a rod-like liquid 
crystal (4-Cyano-4'-n-octyloxybiphenyl, 80CB). 



II. EXPERIMENTAL 

A cuboid (0.4 x 0.43 x 2.65 cm^) was cut from a Vycor 
rod (Corning, code 7930) with a diamond saw. Prior 
to using the Vycor sample and after each imbibition 
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experiment we subjected the host to a cleaning proce- 
dure with hydrogen peroxide and sulphuric acid followed 
by rinsing in Millipore water and drying at 70 ° C in 
vacuum. The cleaned glass was colorless. In agreement 
with former results on the pore structure of samples 
from this Vycor batch , the mean volume porosity, $ 
and mean pore diameter, was determined to $=30% 
and rin =5 nm, resp. 



Water cleaned by a Millipore filtering setup was used 
for the water imbibition experiment. N-hexadccane 
(99.5% purity) was purchased from Aldrich and the 
80CB sample (99% purity) from Synthon Chemicals 
(Germany). The capillary rise experiments were per- 
formed at room temperature, T=22 ° C for water and 
the n-alkane, whereas we have chosen a temperature of 
85 ° C for the imbibition experiment on 80CB. In the 
bulk state, the liquid crystal is at that temperature in 
the isotropic phase. 



The mass uptake of the porous hosts as a function of 
time, m(t), was recorded gravimetrically with an elec- 
tronic balance after gently immersing the monoliths into 
glass vessels containing the bulk liquids. Right upon con- 
tact of the bottom facet of the glass matrix with the 
surface of the bulk liquid, we observed for all samples a 
jump in m{t). It can be attributed to the formation of 
a macroscopic liquid meniscus at the bottom perimeter 
of the matrix, which happened typically within a cou- 
ple of seconds after contact of surface and matrix. Since 
we are interested in the immersion of the liquid into the 
inner pore space, we will present data sets after subtract- 
ing this mass jump in the following. In principle, one has 
also to worry about buoyancy forces acting on the porous 
host and inertia effects in the very beginning of the imbi- 
bition process [H!, However, both contributions are 
negligible in our experiments, a statement supported by 
the data presented below. 
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FIG. 1: Mass uptake of nanoporous silica (Vycor) due to the capillary rise of water (circles), n-hexadecane (triangles), and 
80CB (squares) as function of time, t (upper panel), plotted versus Vt (lower panel). The solid lines in the lower panel 
represent fits of the data sets to the Lucas-Washburn prediction. 



III. RESULTS 



The mass uptake of the porous host by the imbibtion 
of water, 11-C16H32 and 80CB are plotted in Fig. 1 as a 
function of time. For the three investigated hquids we 
find distinct, monotonic increasing mass uptake curves 
and the mass uptake rate, dm(t)/dt is monotonically de- 
creasing with time. Due to the fast imbibition process in 



the case of water, we find additionally to the monotonic 
increase a saturation of the mass uptake at t=27100s and 
ms=m(27100s)=0.136g. At that time the advancing wa- 
ter front reaches the top of the porous monolith. The 
value iRs compares well with the overall mass uptake ex- 
pected for a complete filling of the free pore space of the 
cuboid, 0.136g. 

The time behavior shown in Fig. 1 for the immersion of 
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the three hquids into the porous monoHth is well known 
from capillary rise experiments in tubes and can be de- 
rived by simple phenomenological arguments: Provided 
the liquid wets or partially wets the tube walls (contact 
angle, < 90 ° ), the liquid is sucked into a tube of radius 
r by the capillary pressure, pcap , given by the Laplace 
equation: 

cos 9 

Pcap = 2cr 

r 

The tube geometry calls for a Hagen-Poiseuille type 
flow pattern, where the viscous drag implies that the vol- 
ume flow rate 

dVjt) ^ I pcap 

dt ^ r] l{t) 

and thus the corresponding mass flow rate is given by 

dm{t) 1 Pcap 
dt ^ ^7] l{t) 

, where p is the mass density of the liquid. The variable 
l{t) refers to the height of the liquid column that has al- 
ready been filled up with liquid at time t. A length which 
is at any time proportional to m(t) in a capillary rise ex- 
periment. Thus one arrives at a quite simple relation for 
m{t) and its time derivative: 

, dmlt) 2 Pcap 2^ l^\ 
mit) — ; — p — - ^ p — cosB (1) 
at V V 

In a capillary rise experiment one has also to consider 
gravity forces, however, for the tiny tubes or pores inves- 
tigated here, the gravity pressure is negligible compared 
to Pcap- Relation 1 is solved by 



Since we used for all three experiments the same matrix 
with the same pore microstructure and the same macro- 
scopic extensions we can eliminate the unknown propor- 
tionality factors of relation 1 by a calculation of relative 
imbibition speeds, Wir- The resulting quantities should 
scale according to the ratio of the bulk fiuid parameters, 
provided the continuous-type, macroscopic flow descrip- 
tion as well as the macroscopic fluid parameters are valid 
for the capillary rise phenomenon studied. In order to 
check this, we calculated values of the bulk imbibition 
speeds, vf^^^ by using literature parameters for rj and a 
and contact angles on planar silica. We normalized those 
imbibibition speeds to the value of water: 

From the ratio of the square of the slopes of ■m{t) in 
Fig. 1 (lower panel) divided by p, we extracted the 
imbibition speeds, t;f°"*, in nanoporous silica and, again, 
normalized them by the one of water. As one can 
see in Tab. 1 we get a reasonable good agreement of 
^buik g^j^ij^ v^°'^^ for the liquids investigated. Given the 
uncertainties in the literature values of the bulk fluid 
parameters, which we estimated to be 5-10%, the small 
deviations between wJj""' and vf°'^^a,re well within the 
expected error margins of vjj""^. 

It should also be noted that the decreasing mass up- 
take rate implicates that one probes also a continuous 
range of flow velocities in the pores during a given imbi- 
bition experiment, and thus a continuous range of shear 
rates. Assuming a no-slip boundary condition, one can 
estimate, e.g. for the water experiment, a decrease of 
the shear rate from lO'' 1/s in the very beginning to 1/s 
at the very end of the experiment. Thus, the constant 
slopes of the Lucas- Washburn-law in Fig. 1 as a func- 
tion of time implicate also the absence of any detectable 
shear-thinning or shear-thickening effects for the liquids 
studied. 



m{t) - p^- cosGx/I (2) 

Thus the mass uptake in a capillary rise experiment 
should show a \/t behavior and it is sensitive to the 
fluid parameters (cr, rj, p) and the fluid wall interaction 
(via the contact angle Q). Relation 2 is known as 
Lucas- Washburn law [H, [13] and the square of the 
second prefactor of relation 2 as imbibition speed, fi. 
Porous media, as employed in our experiments, can 
be considered as complex networks of simple tubes. 
Therefore, neglecting details of the structure of our 
porous host, such as the connectivity of the pores, their 
tortuosity and the pore size distribution, we should also 
expect a Vt-behavior for m(t). In agreement with this 
conclusion, the data points of our experiments, plotted 
versus Vi, can be nicely fitted with straight lines - see 
lower panel of Fig. 1. 



IV. SUMMARY 

We have demonstrated that simple capillary rise 
experiments in nanoporous hosts allow one, in principle, 
to study the relative magnitude of fluid parameters of 
liquids conflned in nanopores. In the case of water, 
n-hexadecane and 80CB we observe capillary rise 
velocities that scale as expected from the bulk fluid 
parameters. Thus, we find no hints of a breakdown of the 
continuous- like description of the hydrodynamic flow of 
these liquids in nanoporous silica. This result achieved in 
a capillary flow geometry, and thus in a flow experiment 
with Hagen-Poiseuille type flow geometry, is maybe not 
too surprising, if one recalls the observations reported 
for microscopic thin fllms between mica surfaces in the 
surface force apparatus setup (Couette flow geometry) 
[13, HH Hi- There, deviations from the bulk fluidity 
have been reported for fllm thicknesses on the order of a 
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water 


23 


0.95 


72.25 


1.00 


72.75 





72.75 


1.000 


1.000 


n-Ci6H34 


22 


3.1 


27.15 


0.77 


8.77 


5 


8.77 


0.121 


0.128 


80CB 


85 


9.6 


27.00 


0.98 


2.81 


5 


2.80 


0.039 


0.036 



TABLE I: Bulk fluid parameters and imbibition speeds of the three hquids investigated. The notation is described in the text. 
(The bulk fluid parameters for water were taken from ref. [isl] . the ones of n-Ci6H34 from ref. ^1^] and the parameter of 80CB 
fromref. [13, H El-) 



couple of monolayers and below, only. 

It is understood, that we can miss deviations in the 
absolute transport properties of the spatial confined 
liquids by the type of experiment and analysis presented. 
Our findings allow, however, only two alternative 
conclusions: Either the macroscopic description of the 
nanofluidity in Vycor is valid for all three liquids, or 
it is invalid, but deviates in the same direction and 
probably for a similar reason for all three substances. 
In fact, we assume that the second alternative is the 
correct one. Experiments on the absolute fluid transport 
properties of Vycor glass have been reported in the 
literature [241 , that indicate the existence of one or 
two strongly adsorbed immobile or less mobile layers of 
molecules at the silica pore walls. A finding which is 
reminiscent of the two distinct regimes that one finds 
in sorption isotherms of simple molecules, e.g. argon, 
in Vycor, i.e. an initial layer growth and a subsequent 
formation of pore condensate in the center of the pores 
[1^. Presumably the microscopic boundary layer exists 



also for water, n-Ci6ll34 and 80CB in Vycor and 
decreases or increases the absolute flow rates through 
the pores depending on its mobility and the resulting 
slippage at the pore wall [2, 3, 4, 5]. 

Last but not least, we would like to mention, that 
the imbibition front in Vycor is a fine example of 
an advancing interface in a random environment 
Therefore, imbibition experiments on the capillary rise 
of liquids in this porous host are also interesting in 
order to study kinetic roughening of such interfaces in 
nanostructured environments. 
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